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Abstract. We have used the Stokes radiative transfer code, to model polar- 
ization induced by dust scattering in the polar regions of Active Galactic Nuclei 
(AGN). We discuss the wavelength-dependence of the spectral intensity and po- 
larization over the optical/UV range at different viewing angles for two different 
types of dust: a Galactic dust model, and a dust model inferred from extinction 
properties of AGN. The Stokes code and documentation are freely available at 
http: / /www.stokes-program.info/| 
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Spectropolarimetr y delivers impo rtant results helping to explain the ge- 
Q-f ometry of AGN (e.g. lAntonuccl [2002). The spectropolarimetric properties of 

individual objects are generally complicated and need to be disentangled by 
accurate modeling. The new Monte-Carlo radiative transfer code, Stokes, is 
designed to simulate polarization by scattering in a wide variety of astrophysical 
contexts. The program allows the user to define different emission and scat- 
tering geometries and solve the radiative transfer in 3D. It takes into account 
polarization induced by scattering off free electrons and dust grains. The dust 
composition includes carbonaceous and siliceous grains and the grain size dis- 
tributions parameterized by truncated power laws. The code measures light 
travel times and thus can be used t o model polarization reverberation mapping 
(IShoii. Gaskell. Goosmlmr]|2005l h 



Here we use Stokes to model the polarization expected from scattering in 
dusty, cent rally-illuminated double con es of AGN. We follow the same proce- 
dure as in Goosmann &: Gaskell ( 20071 ). The optical depth of one cone along 



the symmetry axis is ry = 0.3 for the V -band, its half-opening angle equals 
C = 30°. Whilst in lGoosmann fc Gaskelll (|2007l ) we assume a dust composition 



reproducing the extinction in our Galaxy, we here apply a different dust pa- 



rameterization, that inf erred from AGN extinction curves ([Gaskell et al.1 120041 ; 



iGaskell Sz Benker 20071 ). which show a flatter far-UV slope than those for our 



Galaxy. Our AGN dust model implemented here contains 85% "astronomical 
silicate" and 15% graphite. The distribution, n(a), of grain radii, a, obeys to 
the power-law n(a) oc a -2 05 . The resulting polarization and flux spectra, nor- 
malized to the central illuminating flux, are shown in Fig. [1] (left). 
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Figure 1. Modeling a dusty double cone of half-opening angle Be = 30° 
(left) and 9c — 45° (right). Top: polarization, P. Bottom: the fraction, 
F/F*, of the central flux, F*, seen at different viewing inclinations, i. 



Comparison to the results in Goosmann fc Gaskelll ( 20071 ) reveals some clear 
differences: The polarization spectra are weaker at shorter wavelengths than for 
the Galactic dust model. Around 4000 A, the polarization curves drop down 
to nearly zero and then they increase again steeply to longer wavelengths. For 
Galactic dust the polarization spectra show a more gradual increase over the 
whole wavelength range considered. Because of the larger average grain size, 
the scattering efficiency of the AGN dust is significantly smaller. The flux seen 
beyond the edge of the cone, for i > 6c, is therefore low. 

We find similar trends for a larger half-opening angle of Qq = 45° , as shown 
in Fig. Q] (right). However, the obtained polarization percentages are lower, than 
for the narrower cone, because the observer detects the integration of a broader 
range of polarization vectors. The scattered flux at i > 45°, on the other hand, 
is higher due to higher coverage of the central source. 
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